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290 A METHOD OF SOLVING NUMEKICAL EQUATIONS 

similarly x= —2, x= 4 are vertical tangents. The equation\(2X— jw)-4-ju.(\ — 2/x) = 
gives X : ju = 1:1 and 1 : — 1, the directions of the principal axes ; the equation 
p(2X - fi) - X(X - 2ix) = gives X : /* = 2 + a/3 : 1 and 2 - V3 : 1, the directions 
of the asymptotes. 

(3) ' x 2 + 2xy + y 2 - 2x - Qy + 1 = 0. 

df/dx = 2(x -f- y — 1) and df/dy = 2(x + 2/ — 3) have the same direction 
1 : — 1. — df/dx — df/dy = gives x + y = 2, the axis of the curve. Elimina- 
tion of x between / = and df/dx — gives y = 0, the horizontal tangent. 
Similarly x = 2 is the vertical tangent; these tangents cross on the directrix. 
The focus may now be located. 

(4) x 3 + 2tf - Zxy 2 - <ox - Qy = 0. 

df/dx = 3 (a; 2 — y 1 — 2) and df/dy = Q(y 2 — xy — 1). Lines in the two direc- 
tions 1 : 1 and 2 : — 1 cut the cubic in two points only, df/dx + df/dy = 
gives (x— y-\- 2) (x— y— 2) = 0, both asymptotes of the cubic; 2(df/dx) — df/dy— 
gives x 2 -{• xy — 2y 2 = 1, a hyperbola whose center is at the origin and whose 
asymptotes are x — y = and x + 2y = 0, the latter being also an asymptote 
of the cubic. This hyperbola bisects all chords having the direction 2 : — 1. 



A METHOD OP SOLVING NUMERICAL EQUATIONS. 

By S. A. COREY, Hiteman, Iowa. 

The following development of the roots of an equation by Maclaurin's formula 
applies to both algebraic and transcendental equations, and gives all the roots 
approximately whenever the conditions involved in the development can be 
fully complied with. 

Let f(r) = be an equation to be solved, and let a be an approximation to a 
root r of the equation. 

Let us suppose that f(r) is single-valued and analytic in a circle about a as 
a center and including r. Then 

/to = f(a) +f'{a){r - a ) + f -^-(r-a) 2 + f -^ (r - a) 3 + • • •. 

Hence, by the usual formulas for the reversion of series, 1 putting z = f{r) — f(a), 
we get 

r ~ a - \dz ) a Z + \dz 2 ) Q 2l + \dtf )o3\ + "'\dz n ) n\ + '" 

z 2 z 3 d {dr n ~ l \ z n 

where A k = d k f(a)/da k and A = /(a). -In practice it is necessary that the devel- 

1 See, for instance, Goursat, A Course in Mathematical Analysis, Vol. 1, §§ 189-190 (first 
edition). 
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opment in (1) should be rapidly convergent for the number of terms employed, and 
the conditions of convergence for the development make it possible to determine 
the degree of accuracy attained by using a given number of terms. 

It will be observed that of all the derivatives Ak in (1), Ax is the only one which 
has or can have a negative exponent, and hence the only one which must never 
have a zero modulus for any value of a as it approaches r. It will also be observed 
that in order that (1) may be rapidly convergent it is essential that modulus z 
should be as small as possible, and that modulus Ai should be as large as possible. 

Employing only the terms of (1) preceding the term containing z 4 and putting 
2 = f( r ) ~ f( a ) = ~ A, since f(r) = for a root r and A = f(a), we get the 
practical working formula: 

r = a - AAr 1 - U 2 A 2 Ax~ s - W^Mi* + WAiAr* • • • , (2) 

which is quite well adapted to logarithmic computation for both real and imagi- 
nary values of the roots. 

Should the degree of accuracy attained by using (2) not be sufficiently great 
when the first assumed value of a is employed, greater accuracy may be obtained 
by substituting the new value, and so on, until the required degree of accuracy is 
attained. 

To more clearly indicate the use of the method the following examples will be 
useful. 

Let f(x) = x 3 + 2x — sin x — 15 = be an equation, a real positive root of 
which, accurate to five decimal places, is to be found. 

We know from Sturm's theorem, or can learn by trial, graphic methods or 
otherwise, that a real root lies between 2 and 3. Therefore assuming that a = 2, 
we get, 

A = 2 3 + 2 • 2 - sin 2 - 15 = - 3.9093 
A x = 3 • 2 2 + 2 - cos 2 = + 14.4161 

Ai = 3 • 2 • 2 + sin 2 = + 12.9093 

A s = 3 • 2 + cos 2 = + 5.5839 

Substituting in (2), we get x = 2.244. Again letting a = 2.244, finding new 
values of A, Ax, Ai, An, and substituting in (2), we get x = 2.243666, the ac- 
curacy of which is determined not by the convergence of (2) but by the values of 
the sine and cosine as given in the ordinary 6-place tables. 

As a second example let us take the quartic equation 

f(x) = x i - 3a? + 75x - 10,000 = 0, 

which Merriman gives on page 34 of his Solution of Equations, 4th ed., to illus- 
trate Lambert's method. Then taking, as he does, the approximate values of 
the roots to be the four roots of unity each multiplied by 10, we get 
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a = + 10 

A = + 450 
4i = +4,015 
Ai= +1,194 
A s = + 240 



- 10 

-1,050 
-3,865 
+ 1,194 

- 240 



+ 10£ 

+ 300+ 750i 
+ 75-4,060i 
-1,206 

+ 240i 



-lOi 



Ai=+ 24 + 24 + 24 

Substituting in the terms of (2) 
a= +10.00000 -10.00000 



-AAr 1 = - 0.11208 
-iA i A 2 Ax- s = - 0.00186 
-|A 3 ^ 2 2 ^,- 6 = - 0.00006 
+iA s A 3 A l ~ i =+ 0.00001 



= + 9.88601 



- 0.27167 
+ 0.01140 

- 0.00095 
+ 0.00021 

-10.26101 



+ 



300- 750i 
+ 75+4,060i 
-1,206 

- 240i 
+ 24 



Modulus 

a/652,500 

•Vl6,489,225 

1,206 

240 

24 



angle 

68° 11' 55" 

- (88° 56' 30") 

180° 0' 0" 

90° 0' 0" 

0° 0' 0" 



+ 10i 

+0.18330 -0.07728J 

+0.00428+0.00402i 

-0.00012 -0.00033; 

0.00000 • 0.00000 

+0.18746+9.92641; 



-10; 

+0.18330+0.07728; 
+0.00428-0.00402; 

-0.00012+0.00033; 

0.00000 0.00000 
+0.18746-9.92641; 



Should still greater accuracy be required replace the above values of a by the 
values of x just found, and so on until the required degree of accuracy is attained. 



A FORMULA FOR THE SUM OF A CERTAIN TYPE OF INFINITE 

POWER SERIES. 1 

By ELBERT H. CLARKE, Purdue University. 

Introduction. 

The problem to be considered is that of finding a definite, finite formula 
which will give the sum of any convergent infinite series whose terms are such 
that their numerators form an arithmetical progression of any order 2 and whose 
denominators form a geometric progression. 

Since the nth term of an arithmetic progression of the &th order may be 
reduced to the form 

a„ = bon k + hn 10 ' 1 +•■•+&*, 

our problem is to evaluate the expression 

T _ ^b n k + hn*- 1 + ■■■+h 
n=i ar n 



in which the b's are independent of n. But this may be written 
T 






' Y* — - -i- — V - 



a M =i r" 



1 The author wishes to acknowledge criticisms and suggestions from Professors R. D. Car- 
michael and A. C. Lunn. 

2 See Text Booh of Algebra, Chkystall, Vol. 1, page 484. 



